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Abstract. Using the Uves echelle spectrograph at the ESO VLT, we observed the absorption line spectrum of the QSO 0103- 
260 in the Fors Deep Field. In addition to the expected hya forest lines, we detected at least 16 metal absorption systems 
with highly different ionization levels in the observed spectral range. The redshifts of the metal absorption systems are strongly 
correlated with the redshift distribution of the high-z galaxies in the Fors Deep Field and with the strength (but not the number 
density) of the Lyo- forest lines. Both the metal systems and the galaxies show clustering at least up to the QSO emission line 
redshift of 3.365, but only few of these galaxy accumulations seem to form bound systems. 
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1. Introduction 

Observations of high-redshift galaxies and of the absorption 
spectra of high-redshift QSOs are among the most important 
sources of information on the first few Gyrs of the cosmic evo- 
lution. Moreover, since model calculations for the formation of 
galaxies predict that galaxies started to form by gravitational 
contraction of gas concentrated in the deep potential wells of 
the dark matter, at high redshift the intergalactic gas density 
and the galaxy density are expected to be closely correlated. 
The details of this correlation are expected to yield impor- 
tant information on the formation history of galaxies and on 
the interaction of young ga laxies with thei r envir onment. As 
demonstrated recently by Adel berger et alJ d2003l) . the winds 
and the UV radiation of young galaxies can ionize the sur- 
rounding gas and can thus produce gaps in the Lyman forest 
line distribution in the spectra of QSOs. Since atomic nuclei 
heavier than those of Lithium cannot have been formed dur- 
ing the Big Bang, the metal absorption features of the IGM 
contain valuable information on the population conten t and 
star formation history of the early galaxies (see e.g. Songaila 
l200ll IPettini et al]l2002l) . Studies of IGM absorption lines in 
the spectra of QSOs and of high-redshift galaxies along the 
same line-of-sight promise useful information on the questions 
mentioned above. Investigations of this type have been ca rried 
out in the south ern Hubble Deep Field (HDF-S ) bv lTresse et alJ 
dl999l) and by ICristiani & D'Odoricol J2000t). and in the di - 
rection of 8 high-redshift QSOs by Adelberg eretaD |2003). 



In the present paper we describe a similar investigatio n us- 
ing data from the Fors Deep Fie ld (FDF, cf. Appen zeller et alJ 
2000] iHeidt et alJ l200ll l2003t4) . As pointed out in the pa- 
pers quoted above, one of the selection criteria for the FDF 
was the presence of a bright high-redshift QSO (Q 0103- 
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* Based on observations obtained with Uves at the ESO Very Large 
Telescope, Paranal, Chile (proposal No. 66.A-0133). 



260, cf. Iwarren et al.lll99lt iDietrich et ai]ll999t) which was 
included specifically for the purpose of studying such cor- 
relations. Although in the course of the survey 7 additional 
QSOs were detected in the FDF, Q 0103-260 is still the only 
z > 2 quasar which is bright enough for high-resolution spec- 
troscopy. 

The main emphasis of the present investigation is on the 
correlations observed between the galaxy redshift distribu- 
tion and the redshifts of the metal line systems of Q 0103- 
260. However, some information on the Hi absorption towards 
Q 0103-260 is also presented. 



2. Observations and data reduction 

The data presented here are based on observations carried out 
in service observing mode in November and December 2000 
with Uves, the Ultraviolet and Visual Echelle Spectrograph at 
the Nasmyth platform B of ESO's VLT UT2 (Kueyen) on Cerro 
Paranal, Chile. For all observations the red arm of Uves, with a 
central wavelength of 5200 A and a slit width of 1" (resulting 
in a usable spectral range of 4166 - 5164 A and 5230-6212 A 
and a (measured FWHM) spectral resolution of almost exactly 
40 000), was used. The wavelength gap near the central wave- 
length is due to the gap between the CCDs of the Uves red-arm 
detector mosaic. Further technical information on Uves can be 
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Table 1. Observing dates and conditions 



Date 


Start (UT) 


Airm. 


Seeing 


Remarks 


Nov. 17 


02:17:38 


1.00 


07 5-0? 8 




Nov. 17 


03:10:22 


1.03 


075-078 




Nov. 18 


03:19:15 


1.05 


0/4-0.' 8 




Nov. 18 


04:12:13 


1.15 


0/4-0/8 




Nov. 19 


03:19:53 


1.06 


0/6-0/8 




Nov. 19 


04:12:44 


1.17 


0/6-0/8 




Dec. 14 


00:40:34 


1.00 


0/6-173 


moonlight 


Dec. 16 


02:46:52 


1.23 


0/7-0/9 





found in the Uves User's Manual JP'Odorico et al.ll2002l) on 
the ESO Web page. 

The total integration time of about 6.9 hours was split into 
8 individual exposures of 3100 s each. During most nights 
the observing conditions were photometric. One exposure was 
somewhat degraded by moonlight contamination. Basic infor- 
mation on the individual frames is listed in Tabled All origi- 
nal spectroscopic frames are available from the ESO archives 
in Garching. 

In order to obtain some information on the spectral re- 
gion corresponding to the Uves CCD gap, we supplemented 
our high resolution Uves observations with a medium resolu- 
tion (measured FWHM resolution 2 800) spectrum covering the 
spectral range 4530 - 5820 A using Fors2 with Grism 1400V 
at the ESO VET UT2. The S/N per resolution element of the 
medium resolution spectrum is > 100. 

Since the Uves observations were carried out in service 
mode using standard settings, the basic spectroscopic reduc- 
tion procedures (bias correction, background subtraction, cos- 
mic ray removal, flat-field correction, order extraction, sky sub- 
traction, wavelength calibration, merging of the orders and re- 
sampling to a common wavelength scale) were carried out us- 
ing ESO's Uves pipeline reduction procedures. Flux standard 
star observations (with a 10" slit) were obtained once during 
each of the observing nights and were used for a flux calibra- 
tion of the reduced data. 

The identification and subsequent an alysis o f the a bsorp- 
tion lines was based on the catalog of iMortonl (119921) . who 
lists below 2000 A vacuum wavelengths. Moreover, the line 
fitting programme VP FIT was used which assumes vacuum 
wavelengths as input. Therefore, in addition to the heliocen- 
tric radial velocity correction, a transformation from air wave- 
lengths to vacuum wavelengths was applied to all spectra. This 
conversion was ca rried out following th e official IAU rules 
JOosterhofJ i ll 9571) . based on lEdlerJ i fl953k 

The final spectrum resulting from the co-addition of the 
individual reduced spectra and representing a total integration 
time of 24800 s is displayed in Fig.^ It nas a continuum S/N 
per resolution element of about 35 - 50 in the wavelength range 
4200 - 5160 A and about 50 - 65 in the range 5230 - 6210 A. 
The mean error of the wavelength calibration (as derived from 
the dispersion solution) is about 10 mA in the blue section and 
about 7 mA in the red section of the spectrum. 



of theFDF QSO 0103-260 

250 | , , 1 1 , 1 , , 




4500 5000 5500 5000 

Wavelength (A) 



Fig. 1. The complete Uves spectrum of Q 0103-260. The gap 
extending from 5164 A to 5230 A resulted from the space be- 
tween the two CCDs of the Uves red arm CCD mosaic. 



Since this investigation deals with the measurement of ab- 
sorption features, some effort was devoted to determine a suit- 
able and meaningful continuum level and to assess the accu- 
racy of this function. As all absorption lines are very narrow 
compared to QSO emission features, we were not interested in 
the true (non-thermal and thermal) continuum emission of the 
quasar. Instead, we tried to derive the quasi-continuum formed 
by the superposition of the true continua and the broad emis- 
sion lines and blends of the QSO. As this quasi-continuum is 
difficult to predict from models, it was determined empirically 
from local 5th-order polynomial fits through apparently line- 
free sections of the observed quasi-continuum. This procedure 
is relatively safe in the red spectral region, even at the position 
of broad emission lines (except at the peak of the Lycr emission 
line, which was omitted from the following analysis). However, 
due to the large number of Lya absorbers, the assumed quasi- 
continuum blue-ward of the Lya emission line peak depends 
somewhat on the subjective judgement of the reality of appar- 
ent continuum islands between the Ly forest lines. In order 
to estimate the uncertainties possibly introduced by this sub- 
jective judgement, we repeated the procedure of selecting the 
fulcrums and the subsequent interpolation 10 times. At posi- 
tions where the quasi-continuum level appeared uncertain, we 
used, on purpose, different assumptions in the different runs. 
The 10 versions of the quasi-continuum were then averaged, 
and the variance between the different versions was used to 
estimate the accuracy of the continuum derivation. We found 
differences of the local continuum level derived in this way of 
up to 7% in the region of the Ly/3, Ovi and Lyy emission lines. 
In most parts of the spectrum, however, the deviations were a 
few percent at most. Therefore, the uncertainties of the contin- 
uum level, while present, should not have influenced the sta- 
tistical conclusions described in this paper. We also note that 
we checked for, but did not find any evidence for, a possi- 
ble under-cor rected zero-order contamination, as described by 
Rauch ( 1998), which could have influenced our measurements. 
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3. Line profile fits 

Although the present study concerns mainly the investigation 
of the metal lines in the spectrum of Q 0103-260, it was im- 
portant to first get some information on the Lya absorption 
lines which form the great majority o f the absorption features in 
the observed spectrum. According to lCarswell et alJ ll984). the 
profiles of typical Lya lines in high-resolution (FWHM < 25 
km s _1 ) QSO spectra can normally be approximated by Voigt 
profiles. Each line can hence be characterized by its Doppler 

parameter b = ^IkTnT^ + £> 2 urb , the column density /Vhi, and 
the redshift z. 

The line profile 1(A) for each line is then given by 

I(A) = I xe- TW (1) 
where the optical depth t(A) is 

t(A) = x V(a, u) (2) 

b V21n2 

and V(a, u) is the Voigt profile function. 

A powerful and convenient fit program for Ly forest lines 
is the code VPFIT which is based on the above as sumptions 
and which has been developed by ICarswell et al.l dl984l) . It 
is available from (www.ast.cam.ac.uk/~rfc/vpfit.html). This pro- 
gram, which uses /Vhi, b, and z as fit parameters and minimizes 
a suitable ^ 2 function of these parameters, was used throughout 
the present study. 

The main application of VPFIT was the identification and 
fit of the Hi absorption lines in the region blue-ward of the 
Lya emission. As a first approximation, all these lines were 
assumed to arise from hydrogen. Hence, the identified metal 
lines were also included in the fit. Since VPFIT works best with 
small data "chunks", the spectrum was divided into 16 sections 
of roughly 40 A overlapping by up to 5 A. Each section was 
fitted separately. 

Normally, following Ki m et all J200ll) . more lines should 
be added to an existing fit if x\ ^ 1-3. In some cases where 
regions contained either a saturated line or a multitude of very 
weak lines with column densities /Vhi < 10 125 cmT 2 ,xl na d to 
remain greater than the 1 .3 threshold. (Fixing this threshold ob- 
viously is to some extent arbitrary, and increasing the number 
of components normally improves the fit. There is, however, no 
guarantee that a better fit, i.e. a smaller^ 2 , is always physically 
"more correct"). 

Tests with different continuum normalizations and with fil- 
ters applied to the spectrum before the fitting procedure showed 
that the fit results for lines with ~ 13.2 < log /Vhi < 15.5 are 
very robust against such differences in the data handling. For 
very weak lines (log /Vhi < 12.8) and for saturated lines, the fit 
results were found to depend on details of the data reduction. 
Thus, the results for such lines should be taken with care and 
anybody wishing to make use of our results on very weak or 
saturated lines should repeat the data reduction. 

An example of the results of the line fitting procedure with 
VPFIT is given in Fig.|2] 




4705 4710 4715 4720 



Wavelength [A] 

Fig. 2. Example of the spectral fits obtained with VPFIT. 
The normalised observed spectrum and the model spectrum 
(smooth line) are plotted for a 20 A section of the Lyman for- 
est. Also plotted is the difference (for clarity downshifted by 
0.25) between the fit and the data. The narrow lines at ~ 4703 - 
4705 A are metal lines ( Sirvil394 at z ~ 2.375 and Fen/12383 
at z ~ 0.974). The tick marks above indicate the positions of 
the line centroids of the lines found by the fit procedure. For 
the saturated region near 4708 A (where the fit is ambiguous, 
see text) assuming a blend of three lines gave the best result. 



4. The Lye forest lines 

In order to verify that the intergalactic medium in the direc- 
tion of the FDF has normal properties, we compared our new 
data with other Lyman forest studies of QSOs at similar red- 
shifts. For this purpose, we calculated and compared various 
statistical quantities and functions, which are briefly described 
in the following paragraphs. To avoid confusion with the Ly6 
forest, features caused by the strong Ovi lines, and the broad 
blends close to the Lye emission peak, the Lya forest analysis 
was restricted to the wavelength range from 4530 to 5165 A. 
In this wavelength range the fit procedure resulted in the de- 
tection of 423 lines. Of these lines, 49 were removed from the 
sample because their small line widths (b < 15 km s -1 ) made 
an identification with metal absorption more likely, or because 
a too low (/Vhi < 10 12 5 cnT 2 ) or too high (N m > 10 16 cnT 2 ) 
column density made their analysis uncertain. Therefore, only 
374 "reliable" lines were used in the following analysis of the 
Lyman forest. However, since the fraction of discarded lines is 
still small, the concentration on the "reliable" Hi lines should 
not have had any significant influence on the statistical results. 

A plot of the Hi fine fit from 4530 to 5165 A as well 
as a line list derived from this fit is available on the web 
site http://www.astronomy.ohio-state.edu/~frank The correspond- 
ing catalog lists the redshifts z a b s , the column densities /Vhi, 
and the broadening parameters b for the reliably measured Lyor 
absorption lines. 

For the 374 lines selected in this way we investigated the 
following quantities and relations: 
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4.1. Column density distribution function 

The column density distribution function (CDDF), f{Nm), 
measures the number of absorption lines per unit column den- 
sity (in units of cm 2 ) and per unit absorption path (defined as 
x( z ) = 1/2 [(1 + zf - 1.0] for q = 0.0, which was adopted 
here primarily for the the reason of comparability to data sets 
by other authors) as a function of Nm. It is comparable to the 
luminosity function in the study of galaxies. For the CDDF of 
the Lya forest of Q 0103-260 we find that the distribution can 
be well approximated by 



f(N m ) = 7.1 x lOX 1 ' 62 

which is in good agreement with earlier results demonstrat- 
ing that the CDDF can always be approximat ed by a power 
law with an expo n ent near -1 .5 (IHu et all 1 995: Lu et al. 1996; 
Kim et all Il997t iKirkman & Tvtlerl Il997t iKim et alJ 1200 it 
Wevmann et al.ll998HLu et alJl99ll) . 



4.2. Line number density 

For the line number density, defined as 



dN_ 

dz 



dz) 



where dN is the number of lines per redshift range dz, (77)0 
is the local number density of the forest (=a 9. 1 + 0.4), and 
y w 2.19 + 0.27, according to lKim et all J200ll) . we find for 
13.6 < N m < 16.0: 



log(-) = 2.27 ±0.05 
dz 



at a mean z of 2.9 9, which falls exa ctly on the corresponding 
relation derived by Kim et al. (2001) for z > 1.5. 



Table 2. Line list for the three metal line systems with red- 
shifts close to that of the QSO emission lines. N denotes the 
minimum number of velocity components observed in the line 
profile. An asterisk in the last (Re) column points to a com- 
ment in Section 5.2. "bi" in this column indicates that the line 
is significantly blended. 



/lobs 


Ident 




z 


N 


EW 


m.e. 


Re 


(A) 










(A) 


(A) 




z3 


4267.2 


Cm 


977 


3.367 


4 


<.67 




bl 


4506.8 


Ovi 


1032 


3.367 


4 


<2.0 




bl 


4531.7 


Ovi 


1038 


3.367 


4 


<2.0 




bl 


5269.3 


Sim 


1207 


3.367 


4 


<.68 


.038 


* 


5410.4 


Nv 


1239 


3.367 


4 


1.503 


.032 




5427.8 


Nv 


1243 


3.367 


4 


1.222 


.020 




5827.9 


Cn 


1335 


3.367 


4 


0.032 


.015 




6086.6 


Siiv 


1394 


3.367 




<0.05 






6125.9 


Siiv 


1403 


3.367 




<0.05 






z2 


4262.8 


Cm 


977 


3.363 


4 


<2.6 




bl 


4318.5 


Nm 


990 


3.363 


7 


0.329 


.060 




4502.4 


Ovi 


1032 


3.363 


7 


<2.9 




bl 


4527.2 


Ovi 


1038 


3.363 


7 


<3.8 




bl 


5264.1 


Sim 


1207 


3.363 


7 


1.190 


.066 




5405.1 


Nv 


1239 


3.363 


7 


2.072 


.057 




5422.5 


Nv 


1243 


3.363 


7 


1.444 


.061 




5422.6 


Cn 


1335 


3.363 




<.04 






6081.1 


Siiv 


1394 


3.363 


7 


1.280 


.072 




6120.3 


Siiv 


1403 


3.363 


7 


0.906 


.065 




zl 


4256.0 


Cm 


977 


3.356 


4 


0.920 


.029 




4494.9 


Ovi 


1032 


3.356 


5 


<4.0 




bl 


4519.8 


Ovi 


1038 


3.356 


5 


<4.2 




bl 


5255.6 


Sim 


1207 


3.356 


5 


0.460 


.035 




5396.4 


Nv 


1239 


3.356 


5 


0.490 


.042 




5413.8 


Nv 


1243 


3.356 


5 


0.290 


.037 




5813.2 


Cn 


1335 


3.356 




<.03 






6071.2 


Siiv 


1394 


3.356 


5 


0.430 


.028 




6110.5 


Siiv 


1403 


3.356 


5 


0.290 


.030 





4.3. Mean HI opacity 

For the mean or effective Hi opacity r e ff , defined as 
e -r dT =< e -r > ( wnere <> indicates averaging over the wave- 
length) we obtain for z = 2.99 the value 0.36 ± 0. 05, w hich 
is aga in in excellent agreement with the results of IKim et alJ 
(2001) and others. 

4.4. Line counts 

The number of lines in a given z interval as a function of the flux 
threshold F u and the number of lines as a function of the filling 
factor, defined as the fraction of the spectrum occupied by pix- 
els whose normalized flux is smaller than a given flux threshold 
F t , again, were found to a gree with the corresponding functions 
derived in earlier studies (IKim et alJl200lUSchave et al.ll 19991 
OutrarrJl994lKirkrrian & Tvtleill997UHu et alll995UFu et al] 
1996:IPress et al.ll993l). 



4.5. Doppler parameter distribution and 
temperature-density relation 

The distribution of the Doppler parameter b for the Lya 
lines can be well approxim ated by a Hui-Rutledge distribu- 
tion llHui & Rutledgefl 19991) with b H R = 13.8 km s" 1 and b a = 
(23.7 ± 2.8) km s _I . In order to derive information on the tem- 
perature of the IGM, we investigated the relationship between 
the (baryon) density and the temperature T by looking for a 
cut-off in the observed column density vs. Doppler parameter 
distribution. 

The resulting temperatures of ~ 15 000 - 20000 K at the 
mean column density of Nm =1.9 xlO 13 cirT 2 are in good 
agreement with photoionization models. 
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Table 3. Line list for the metal absorption line systems with 
2.1<z<3.2 



dobs 

(A) 


Went 




z 


N 


EW 

(A) 


m.e. 

(A) 


Re 


4300.5 


Ovi 


1038 


3.145 


2 


0.045 


.020 


* 


5134.6 


Nv 


1239 


3.145 


2 


0.047 


.021 


* 


5151.1 


Nv 


1243 


3.145 


2 


0.041 


.020 


* 


5531.2 


Cn 


1335 


3.145 


2 


0.019 


.020 


* 


5776.7 


Siiv 


1394 


3.145 


2 


0.194 


.029 




5814.2 


Siiv 


1403 


3.145 


2 


0.130 


.025 




4271.2 


Ovi 


1038 


3.116 


2 


0.029 


.020 


* 


4966.5 


Sim 


1207 


3.116 


2 


0.347 


.035 


* 


5115.6 


Nv 


1242 


3.116 


2 


0.021 


.012 


* 


5737.3 


Sirv 


1394 


3.116 


2 


0.154 


.023 




5774.4 


Siiv 


1403 


3.116 


2 


0.120 


.022 




4557.4 


Nv 


1239 


2.679 




<0.04 


0.015 


* 


5695.6 


Civ- 


1548 


2.679 


! 


0.105 


.013 


* 


5705.0 


Civ 


1551 


2.679 


1 


0.063 


.013 


* 


5508.6 


Civ 


1548 


2.558 




0.268 


.023 




5517.8 


Civ 


1551 


2.558 


[ 


0.169 


.022 




4958.3 


Siiv 


1394 


2.558 


l 


0.042 


.020 




5507.8 


Civ 


1548 


2.558 


[ 


0.155 


.021 




5517.0 


Civ 


1551 


2.558 




0.090 


.013 




4703.3 


Siiv 


1394 


2.375 


2 


0.430 


.030 




4733.7 


Siiv 


1403 


2.375 


2 


0.385 


.040 




5224.4 


Civ 


1548 


2.375 


2 


0.350 


.068 


* 


4330.6 


Cii 


1334 


2.245 


2 


0.188 


.010 


* 


4522.8 


Siiv 


1394 


2.245 


2 


0.350 


.020 


* 


4552.0 


Siiv 


1403 


2.245 


2 


0.290 


.015 


* 


4954.2 


Sin 


1527 


2.245 


2 


0.068 


.025 


* 


5023.9 


Civ 


1548 


2.245 


2 


? 




bl 


5032.3 


Civ 


1551 


2.245 


2 


? 




bl 


6018.6 


Aim 


1855 


2.245 


2 


0.109 


.010 




6044.8 


Aim 


1863 


2.245 


2 


0.045 


.010 




4329.6 


Cn 


1334 


2.244 


2 


0.042 


.014 




4521.7 


Siiv 


1394 


2.244 


2 


0.112 


.020 




4551.0 


Siiv 


1403 


2.244 


2 


<.ll 




bl 


5022.7 


Civ 


1548 


2.244 


2 


<.19 




bl 


4421.8 


Siiv 


1394 


2.173 


1 


0.133 


.010 




4450.5 


Siiv 


1403 


2.173 


[ 


0.095 


.011 




4911.8 


Civ 


1548 


2.173 




0.360 


.020 




5884.2 


Aim 


1855 


2.173 


j 


0.032 


.015 




5909.6 


Aim 


1863 


2.173 




0.042 


.016 




4910.8 


Civ 


1548 


2.172 




0.180 


.015 


* 


4919.0 


Civ 


1551 


2.172 




0.157 


.016 


* 


4836.2 


Civ 


1548 


2.124 




0.168 


.010 




4844.2 


Civ 


1551 


2.124 




0.122 


.012 




5793.6 


Aim 


1855 


2.124 




0.031 


.019 




5818.8 


Aim 


1863 


2.124 




0.015 


.015 





Table 4. Line list for the two metal absorption line systems near 
z = 0.974 



'lobs 


Went 




z 


N 


EW 


m.e. 


Re 


(A) 










(A) 


(A) 




4688.1 


Fen 


2374 


.9744 


5 


0.260 


.080 


bl 


4704.5 


Fen 


2383 


.9744 


5 


0.399 


.030 


bl 


5107.0 


Fen 


2587 


.9744 


5 


0.240 


.035 


bl 


5133.7 


Fen 


2600 


.9744 


5 


0.438 


.030 


bl 


5521.1 


Mgn 


2796 


.9744 


5 


1.080 


.040 




5535.2 


Mgn 


2804 


.9744 


5 


0.880 


.030 




5632.8 


Mgi 


2853 


.9744 
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From the above comparisons we conclude that the Lya for- 
est of the FDF quasar has parameters which are rather typical 
for this redshift and that the IGM along the line of sight to 
Q 0103-260 has normal properties. 

5. The metal lines 

5.1. Line identifications 

About 10 % of the absorption lines in the Uves spectrum belong 
to UV transitions of common ions of elements heavier than He 
in various ionization states. Following astrophysical traditions, 
all such lines are called "metal lines" in this paper. Except for 
very weak and very strong lines, the metal lines in the Uves 
spectrum are easily recognized on the basi s of their lower line 
width (see e.g. Fig. As suggested bv iRauch et alJ (fl997b) 
we assumed all reasonably strong lines with b < 15 km s _I to 
be metal lines. In order to assign these lines to certain atoms, 
ions, and redshifts, we used the following procedure: First, we 
identified the absorption lines connected with the QSO and its 
immediate environment. From earlier studies of the Fors Deep 
Field it is known that Q 0103-260 is embedded in an accumu- 
lation of bright sta rburst galaxies at nearly the same redshift 
( iHeidt etaljl2003al) . Hence, we expected (and found) strong 
absorption line systems superimposed on the quasar emission 
lines. Because of their association with the quasar emission 
lines and their known redshifts (within 0.01 of the emission 
line redshift) these lines could be readily identified indepen- 
dently of their width. In total we found three different groups 
of lines with z ~ Zqso- In Table [2] these redshift groups are 
labeled (in the order of increasing redshift) as zl, z2, and z3. 

As a next step, we used our identified Lya forest lines to 
look for metal emission at the redshifts of the strongest un- 
blended Hi absorption lines. The rest wavelengths and expected 
line strengths o f the metal lines were taken from the catalog of 
iMortonlll 19921) . For those narrow lines which were still uniden- 
tified after this step, we calculated mutual wavelength ratios 
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which were compared with those expected for common Ly for- 
est metal lines (making use of the fact that wavelength ratios 
are invariant to the redshift). At the high spectral resolution of 
Uves, this method is particularly powerful for identifying lines 
belonging to multiplets. For each line we then calculated the 
redshift. A metal absorption system was assumed to be estab- 
lished when at least 3 lines or blends (in some cases including 
an Hi line) were found to have redshifts not deviating by more 
than 0.0001. Having established a metal system in this way, we 
searched for other expected lines at the same redshift. 

Using this procedure, we were able to identify more than 80 
metal lines or blends belonging to 16 redshift systems. About 
20 (mostly weak) sharp lines remained unidentified. Lines were 
assigned to the same redshift system if their profiles overlapped 
for at least one metal line. With few exceptions each system 
contains several components. In some cases an optimal line 
fit required the assumption of up to 16 components. However, 
as pointed out above, optimized line fits may overestimate the 
number of real components. Therefore, we determined for each 
system also a lower limit of components by counting the pro- 
file components which were directly visible on our spectra. 
The metal lines and blends, for which reliable observed wave- 
lengths or equivalent widths could be measured, are listed in 
Tables |2] to Lines and blends belonging to the same system 
are grouped together. In cases where only one component of a 
doublet is listed, the second component could not be measured 
because of blending. But in all these cases it has been verified 
that the spectrum is consistent with the presence of both com- 
ponents at the expected absorption strength ratio. The individ- 
ual columns of the tables give the mean observed wavelength of 
the line (averaged over the individual velocity components), the 
line identification, the approximate mean redshift z, the lower 
limit of velocity components N, the total observed equivalent 
widths EW (derived by summing up the contributions of all 
individual velocity components), and their mean errors. An as- 
terisk in the last column indicates a comment in the Subsection 
5.2., whereas "bl" denotes blending with other lines (usually 
Hi). For some lines blending made it impossible to derive re- 
liable equivalent widths. We then give as upper limits the EW 
which the line would have if the contribution of other lines were 
negligible. In addition to data for reliably detected lines, some 
upper limits or non-significant EW values have been included 
in the tables if these values provide information on line ratios 
important for the derivation of the ionization parameter or the 
plasma density. For many redshift systems, additional lines are 
visible in the Uves spectrum, but cannot be measured because 
of blending. This applies in particular to wavelengths below 
4500 A where the Lya and Ly/? forests merge. 

5.2. Remarks on individual metal absorption systems 
z= 3.367 (z3): 

At least part of the absorption at the expected position of the 
Sim resonance line seems to be due to Lya. Therefore, only 
an upper limit for the Sim absorption can be given. Since Siiv 
is absent in this system, and since the strength of Ovi and Nv 
suggest a very high ionization parameter, a significant Sim ab- 



sorption is not to be expected, but cannot be ruled out by our 
observations. 

z = 3.145: 

All metal lines of this system are weak. Only the lines of Siiv 
are significant detections. 

z = 3.116: 

Weak Oiv and Nv and strong Sim indicate a low ionization 
state. Cn, however, is not detected. Sim could be overestimated 
due to blending with a weak Hi line. 

z = 2.679: 

Except for the three lines listed, all lines coincide with regions 
of strong and complex blends. 

z = 2.375: 

EW(CiV/U548) has been measured in the Fors spectrum. 
NV/U243 is clearly detected, but not measurable because of too 
severe blending. 

z = 2.245: 

The system consists of more than two (at least 3, possibly 7) 
components. The wavelength and EW values refer to the two 
strongest and best defined components only. The line ratios in- 
dicate a much lower ionization parameter than observed in the 
other systems. 

z = 2.172: 

The only statistically significant metal lines are those of the 
Civ resonance doublet. Siiv is not detected, but the well de- 
fined sharp lines and the excellent agreement with the theoreti- 
cal wavelength ratio makes it very likely that the identification 
of this weak system is correct. 

z = 0.9737: 

The Feii/12383 A line could not be measured accurately be- 
cause of an uncertain continuum level at the corresponding 
wavelength. The strong lines of low ionization stages indicate 
that the ionization parameter is much lower than in the systems 
with higher redshifts. 

5.3. General properties of the observed metal systems 
5.3.1. The metal absorbers with z abs ~z em (QSO) 

As mentioned above and listed in Table|2] there are three strong 
absorption systems with lines of Cm, Nm, Nv, Ovi, Sim, and 
Siiv, at redshifts clo se to the qu asar's emission line redshift (z = 
3.365, according to Noll 2002). Such systems with z a bs ~ z em , 
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which are observed in many quasar spectra, and are often re- 
ferred to as the "intrinsic" absorbers, have been proposed in the 
literature to be either due to matter ejected from the QSOs or 
to originate in galaxies physically associated with the quasars. 
In a few cases evidence for a close association of such systems 
with the central QSO has bee n reporte d, e.g. time-variability of 
the line strength ( Hamann et aljl9 95|) or partial coverage of th e 
background source (Pet itiean et al.ll99 4: Gang uly et alJl 999). 
In our case, the known accumulation of galaxies with abo ut the 
same redshift around Q 0103-260 (e.g. He idt et alJ l2001) sug- 
gests that galaxies associated with the QSO at least contribute. 
On the other hand, of the three galaxies with reliable spectro- 
scopic redshifts z ~ Zqso and closer than 5 arc sec to th e line 
of sight to the QSO only one (object 9003 in lNollll2002l) has a 
redshift coinciding within the error limits with one of the sys- 
tems listed in Table [3] (z3), while the two other galaxies have 
slightly larger (by 0.008-0.013) redshifts, indicating that they 
are located behind the QSO. However, there are additional faint 
galaxies close to the QSO without known spectroscopic red- 
shifts, which could well cause the other z « zqso absorption 
systems. 

For the system z3 we find a slightly (by about 0.002) higher 
redshift than derived for the QSO emission lines. However, as 
this difference is within the errors of the derivation of the emis- 
sion line redshift, this result cannot be regarded as evidence for 
a motion of the absorbing gas towards the QSO. 

In Figs. and 0] the profiles of the three intrinsic absorp- 
tion systems zl to z3 and the corresponding spectral fits are 
presented. For these fits it was assumed that the same veloc- 
ity components with exactly the same individual redshifts were 
present at all lines, while the absolute and relative strength of 
the components was allowed to vary from line to line. The 
wavelengths of the individual components are, if present, in- 
dicated in Figs. [5]and |4] A table of the fit parameters and the 
column densities derived from the fits is also available elec- 
tronically from the web page noted above. It is evident from 
the plotted profiles that each system consists of discrete ab- 
sorptions with different velocities and that the different veloc- 
ity components have different strengths for different lines. As 
demonstrated by the figures and Table|2] all three systems show 
very strong high-ionization lines while significant absorption 
features of low-ionization stages (like the resonance lines of 
Cn, Nn, Sin, and Oi) were not detected. Siiv is fairly strong in 
zl and z2, but not detected at z3. On the other hand, an ab- 
sorption feature is present at the expected position of the Sim 
line of the z3 system. However, the profile of this feature indi- 
cates that at least part (and possibly all) of this feature is due 
to Lyg. A comparison with photoioni zation models (see e.g. 
Berg eron & Stasinskalfl986l Isteidellfl990l ISongaila & Cowid 
1996) indicates that the ionization parameter T (while high for 
all three "intrinsic" systems) increases from zl to z3. Hence, 
it appears likely that the apparent Sim line of z3 is in fact due 
to Hi at a lower redshift and that Siiv (and lower Si ionization 
stages) have been destroyed by strong radiation short-ward of 
the He + break. 

Because of the relatively few significant line strength mea- 
surements our data do not allow to derive reliable quantitative 
values for the densities, ionization parameters and abundances 
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Fig. 3. Line profiles of the blends produced by the metal ab- 
sorption systems zl to z3 with redshifts close to the quasar 
emission line redshift (thin line). Also included is the the model 
spectrum produced by the fit procedure (smooth thick line) and 
the mean error of the observed spectrum (thin line close to the 
zero level). The solid tick marks above the spectra indicate the 
centroids of the different line components assumed for the fit. 
Dashed tick marks represent the centroids of Hi at other red- 
shifts added to obtain adequate fits. The relatively large resid- 
uals near Lyy are due to a detector defect (at 4249 A) and the 
omission of lower-z Lyo- lines in the fit. 

in the observed metal systems. From the EW ratio Nv/Siiv it 
is clear, however, that the lines are formed in a plasma of rel- 
atively low density and a very high ionization parameter, and 
that the ionization parameter is decreasing with decreasing red- 
shift. Such a behavior is to be expected if the redshift is related 
to the distance and if the QSO provides the ionizing radiation. 



5.3.2. The metal absorbers with 2.0 < z < 3.2 

Compared to the "intrinsic" systems all observed metal sys- 
tems at lower redshifts have line ratios indicating lower ioniza- 
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Fig. 4. Line profiles of the blends produced by the metal ab- 
sorption systems zl to z3 with redshifts close to the quasar 
emission line redshift, continued. The Hi blends between the 
three Sim features could not be fitted unambiguously and, 
therefore, were excluded from the fit. 



tion parameters. Most of these systems, though, are dominated 
by high-ionization lines too, with lines from neutral or single 
ionization stages not being detectable. Only one system (z = 
2.245) shows significant absorption by low-ionization lines. 

The high state of ionization of the metal systems listed 
in Table El agrees well with the theore tical predictions of 
iRauch et all dl997al) and lHaehnehl ill 9971) who used hydrody- 
namic computations to calculate the spatial distribution and 
line absorption of low-metallicity intergalactic gas at high red- 
shift. From a detailed comparison of the relative l ine stren gths 
in TableEJwith the calculated column densities of R auch et al] 
dl997ah - we find the high-ionization lines (such as Ovi, Nv, 
Civ) to follo w the pr edictions rather well. The theoretical re- 
sults of Rauc h et alJ (U£97a) also provide a plausible expla- 
nation for the weakness of absorption lines of low ionization 
species such as Cn, which are expected to be present only in the 
rare cases where the line-of-sight passes through a high density 
core of a protogalactic clump of the intergalactic gas. 

The agreement with the predictions of Rauc h et alJ I ill997ah 
is less good for species with intermediate ionization potentials. 



In particular, we find the Siiv and Sim lines listed in TableGJto 
be on average stronger and the resulting Siiv/Civ and Siiv/Nv 
column density ratios about 10 times higher than predicted for 
a metal-poor gas with a sol ar relative compos i tion. A s in the 
case of the measurements of Son gaila & Co wie ( 1996), the dis- 
crepancy is smaller if our data are compared to the calculations 
made assuming the chemical composition of galactic metal- 
poor stars. But even in this case, our Siiv lines are on average 
too strong by about a f actor of 4. 

A comparison with lMisawa et alJ (l2002h indicates that the 
observed number density of Civ systems with EW > 0.15 A 
in the redshift range 2.1 < z < 3.2 along the line-of-sight to 
Q 0103-260 is normal for a QSO at this redshift. Although ad- 
ditional systems with smaller EW values are present in Table|5J 
our data provide no evidence for a high abundance of weak 
metal systems. 

5.3.3. The Mg/Feu systems at z « 0.974 

While there are several sharp lines in the red part of the spec- 
trum which are suspected to be due to low-ionization lines at 
moderate redshifts, only two metal systems with redshifts < 2 
could be reliably identified. The profiles of the strongest lines 
of these two systems with almost identical redshifts (z ~ 0.974) 
are presented in Fig. [5] Both systems consist of various com- 
ponents and show a complex structure. The plasma produc- 
ing these lines has a relatively low ionization level, compara- 
ble to that of most of the interstellar matter in our galaxy and 
the gas observed in low- z metal absorption systems (see e.g. 
iBergeron & Boissell99lh . 

6. Correlations 

Because of the relatively small observed Lya redshift range 
and the presence of the superimposed metal lines, the distri- 
bution of the Lya absorption-line redshifts provides only lim- 
ited information on the cosmic structure along the line-of-sight 
to the observed quasar. In particular, we found for the red- 
shift range 2.7 < z < 3.3, where we have the most reliable 
data on the Lya absorption, a number density distribution of 
the Lya lines which is almost fully consistent with a random 
distribution of the absorbers along the LOS. A clustering sig- 
nal could be detected only for the two-point correlation func- 
tion in the lowest velocity ran ge (Av ^50 km s~') and in 
the step optical depth function dCen et alJll99a) for velocity 
separations ~ 5-10 km s -1 . The clustering signal is higher 
for the systems with higher column density. This result is in 
good agreement with earlier reports on t he clustering of Lyman 
forest lines on small scales (see e.g. IWeb ^[l987j: |Huetalj 
|l995 [lLu et alJl996tlchernomordikll995tlCristiani et alll995t 
Il997l) . However, while little structure is evident in the number 
density distribution of the Lya absorbers, structure is indicated 
in the total Hi absorption as a function of redshift (cf. Fig.0. 
More specifically, we find an absence of high column density 
absorption and a lower-than-average total Lya absorption in the 
redshift interval 2.8 < z < 3.1. 

More interesting information on the distribution and the 
properties of the absorbing gas along the LOS to Q 0103-260 is 
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Fig. 5. The two Mg/Fen absorption systems at z » 0.974. Also 
plotted is the model spectrum (thick smooth line). The tick 
marks again indicate the positions of the centroids of the lines 
assumed for the fit. The two strong absorption lines at 4703.1 
A and 4703.5 A in bottom panel are due to Siiv/11394 at z = 
2.375. 



provided by the redshift distribution of the metal-line systems 
and by a comparison of this distribution with the galaxy red- 
shift distribution in the FDF. From Tables |3] to @] it is obvious 
that the metal absorption systems also tend to cluster on small 
scales. Moreover, no metal absorption system is detected in the 
redshift range 2.8 < z < 3.1 where the Lya absorption is weak. 
A comparison between the redshifts of the individual metal ab- 
sorption systems and those of the Lya absorption lines shows 
that (as expected) detected metal line systems always corre- 
spond to strong Lya absorption. On the other hand, there are 
strong Lya absorption features without detectable correspond- 
ing metal absorption. A notable example is the Lya absorption 
at z — 2.758. The total Lya absorption of this system is simi- 
lar to or slightly stronger than that for the two systems at z = 
2.558, but no metal absorption is detectable. The upper limits 
for the Civ and Siiv absorption are by factors of at least 0.05 
to 0. 1 lower than the corresponding metal line strengths in the 
z = 2.558 systems. We regard these large differences as clear 



evidence for strong local variations of the metal enrichment at 
this redshift. 

The existence of photometric and spectroscopic redshifts 
for a significant fraction of the FDF galaxies allows us to com- 
pare the redshift distributions of the galaxies with those of the 
absorption line systems discussed here. At presen t photometric 
redsh i fts are known fo r about 6500 FDF galaxies ( B ender et alJ 
120011: ICabasch etalJbOOfl while so far (accurate) spectro- 
scopic redshifts are available for only about 270 galaxies JnoI3 
120021: iNoll et alJEooll) . Some properties of the QSO absorp- 
tion line distribution (such as the lower-than-average density 
of objects at 2.8 < z < 3.1 and an over-density near the red- 
shift of Q 0103-260) are also indicated in the distribution of the 
galaxy's photometric redshifts. However, the accuracy of pho- 
tometric redshifts is much lower than that of the Uves absorp- 
tion line redshifts and generally too low for a meaningful de- 
tailed comparison. The available spectroscopic galaxy redshifts 
(m.e. about 0.002) is adequate for the present purpose, but the 
number of available redshifts is uncomfortably small for sta- 
tistical purposes. Fortunately, the FDF spectroscopic observing 
program is biassed towards high redshift objects. Therefore, for 
the redshift range 2.0 < z < 3.4, corresponding to (most of) the 
metal line redshifts in this investigation, the present spectro- 
scopic catalog is » 50 % complete for a limiting I-Magnitude 
of 24.0. At least for this z range the present spectroscopic red- 
shift catalog can be regarded as representative sample which 
can be used for the present statistical purposes. 

Fig. [6] presents a histogram of the z distribution of the 
presently available accurate spectroscopic redshifts 2.0 < z < 
3.7 in the FDF. Also indicated in this figure are the redshifts 
of the metal absorption systems of Q 0103-260 (broken ver- 
tical lines) in this range. In order to avoid meaningless small 
numbers and in view of the expected peculiar velocities of the 
galaxies, we used histogram bins of Az = 0.01. If the redshifts 
of the 77 galaxies included in Fig.[6]were distributed randomly 
we would have to expect on average 0.45 objects per bin and 
a Poisson distribution of the actual number of objects per bin. 
While the population of most bins is roughly consistent with 
this assumption, various bins show highly significant (corre- 
sponding to 3 to 6 <x) deviations. In particular, the clustering 
of galaxy redshifts at the z values 2.245, 2.34, 2.37, 2.48, 3.14 
and 3.38 cannot be explained by statistical fluctuations of a ran- 
dom distribution. Also statistically significant is the total ab- 
sence of objects with redshifts 2.78 < z < 3.03. Such strong 
deviations from a random distribution of the galaxy redshifts, 
which have been found and discussed before in the context of 
other deep fields (see e.g. lCohen et alJ 1 996l lSteidel etaiT l 998 ; 
Giavalisco 2002) are obviously the result of the cosmic large- 
scale structure. In most cases the galaxies included in a single 
redshift bin are distributed over the whole FOV of the FDF. In 
general, the velocity dispersion in the redshift clumps relative 
to the (minimum) observed size of these structures are too large 
to form bound systems, although in some bins (such as in the 
strong peak near z = 2.37) many objects have the same redshift 
within the error limits, making it likely that these accumula- 
tions are or are becoming bound and are precursors of galaxy 
clusters. 
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As shown by Fig. [6] 13 of the 16 metal system redshifts 
coincide with populated galaxy redshift bins. In the absence 
of a correlation (i.e. assuming random distributions) less than 
4 such coincidences are expected to occur. Hence, Fig. ^pro- 
vides not only strong evidence for the presence of large-scale 
structure in the galaxy redshift distribution, but also for a strong 
correlation of the galaxy distribution and the metal absorbers 
which obviously trace the same large-scale structure. The fact 
that no galaxies have been found at the redshifts of the metal 
absorptions system with redshift 3.116 and the two adjacent 
systems near z = 2.558 can readily be explained by the incom- 
pleteness of our spectroscopic survey. More interesting is the 
fact that for the conspicuous clustering of galaxy redshifts near 
Z = 2.34, in spite of a thorough search at the expected wave- 
lengths, no trace of corresponding metal absorption systems 
could be detected. Whether Lya absorption is present at this 
redshift, cannot be tested with our data since the correspond- 
ing wavelength is outside our Uves spectrum. Interestingly, the 
galaxies belonging to this clustering differ from the other red- 
shift groups by a pronounced clustering tendency on the sky, 
with the center of their distribution about 2 arcmin off the LOS 
to the QSO. 

For the small spectral range, where simultaneous informa- 
tion on the Lya lines, the metal lines, and the galaxy redshift 
distribution is available, we have included the data for all three 
components in Fig.0 For this purpose, the Uves spectrum has 
been smoothed to a redshift resolution of 0.005 to be compa- 
rable with the resolution of the galaxy redshifts. Moreover, the 
Uves spectrum has been normalized to the likely continuum 
level in the unsmoothed spectrum. The histogram of the galaxy 
redshift distribution has been scaled by a factor 0. 1 . The figure 
clearly shows the lower-than-average Lya absorption and the 
absence of metal line systems and galaxy redshifts in the z in- 
terval 2.80 < z < 3.05 (corresponding to a comoving distance 
of about 300 Mpc). It is also interesting to note that the strong 
Lyo- absorption near z = 2.76 (which does not show correspond- 
ing metal absorption), while not exactly coinciding with the 
redshift of a galaxy spectrum, nevertheless falls into a redshift 
region where galaxies are observed. 

We also checked our data for evidence of the weakening 
of Lya absorption at the redshift of st arburst galaxies close to 
the LOS reported by Adelbergeret al. (2003). Although some 
of our foreground galaxies are close enough to potentially pro- 
duce such an effect, our numbers are too small to provide reli- 
able conclusions on this question. 

As noted above, only two metal systems with (almost iden- 
tical) redshifts < 2.0 could be reliably identified in our Uves 
spectrum. Fig.[8]shows the distribution of the presently avail- 
able spectroscopic redshifts and the position of the two (on 
this scale coinciding) metal system redshifts. When interpret- 
ing Fig.|8| one has to keep in mind that (in contrast to Fig.[6ji 
the overall redshift distribution in Fig.[8]is strongly affected by 
an observational (negative) bias affecting galaxies in the red- 
shift range 1.2 - 2.0 (cf. Noll 2002). On the other hand, the 
small scale structure of the redshift distribution in Fig.[8]should 
be representative. As illustrated in Fig. [8] the redshifts of the 
two detected low-z systems again coincide with a clustering of 
galaxy redshifts. However, the plot also shows that there are 
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Fig. 6. The redshift distribution of FDF galaxies with accurate 
spectroscopic redshifts in the range 2.0 < z < 3.7. The bro- 
ken vertical lines indicate the redshift positions of the metal 
absorption systems falling into this range. Note that the close 
groups of metal systems near the redshifts 2.172, 2.244, and 
2.558 (cf. Tabled are not resolved on this scale and appear as 
single systems in this figure. 

many other such clusterings at low redshifts for which no ab- 
sorption lines were detected. A reason that just the systems at 
Z = 0.974 are observed could be the fact that some galaxies 
with this redshift are exceptionally close to the line of sight to 
Q 0103-260. One galaxy falling i nto the correspond ing redshift 
bin (No. 4442 in the catalog of iHeidt etai1l2003bl) has an an- 
gular distance of only about 11". While it still appears unlikely 
that gas of this galaxy is responsible for the observed metal 
absorption, there are several galaxies with similar photometric 
redshifts (but no spectroscopic data) even closer to the LOS 
to the QSO, which may be associated with FDF 4442. These 
galaxies appear to be relatively normal late-type galaxies which 
could well produce the observed redshifted interstellar absorp- 
tion features. 

7. Conclusions 

Our analysis of the absorption lines of the quasar Q 0103-260 
and a comparison with spectroscopic redshifts of FDF galax- 
ies confirm that the absorbers' redshift distribution traces well 
a common large scale structure. The fact that the number den- 
sity of the Lya lines shows no significant deviation from a ran- 
dom distribution in the observed redshift range, while the total 
Hi absorption shows variations, seems to indicate that the Ly 
forest clouds in the region investigated show systematic local 
variations either of their density or of their ionization level, cor- 
related with the galaxy density. The presence of Lya clouds of 
high neutral Hi column densities in the FDF region is found to 
be highly correlated with strong metal line absorption and the 
presence of starburst galaxies at the same redshift. This obvi- 
ously suggests that hydrogen clouds with high neutral column 
densities are (or have been) places of early vigorous star for- 
mation and the enrichment of the corresponding clouds with 
heavy elements. In agreement with other authors, we find large 
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Fig. 7. The normalized flux in the QSO spectrum showing the 
strength of the Lya absorption as a function of redshift. Also 
given are the redshift positions of the metal absorption sys- 
tems in this redshift range (vertical broken lines) and the his- 
togram of the redshift distribution of the FDF galaxies (scaled 
by a factor of 0.1). For this figure the Uves spectrum has been 
smoothed to a resolution of Az = 0.005. It has been normalized 
to the continuum level of the unsmoothed spectrum. Note that 
the strong broad absorption features at z = 2.70 and 2.72 are in 
fact due to Ovi absorption at z = 3.36. The Lya absorption at 
2.68 is contaminated by Ly/3 at z = 3.36 and Lyy of this redshift 
contributes to the feature at 2.49. Most other strong features are 
due to Lya. 
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Fig. 8. Distribution of spectroscopic galaxy redshifts in the 
FDF for 0.5 < z < 2.0. The relatively low frequency of objects 
between 1.0 and 2.0 is due to an observational selection effect. 
The vertical broken line indicates the position of the z = 0.974 
Fen/Mgi/Mgn absorption system described in the text. 



variations in the metal line strengths and ionization parameters 
of individual Lya absorbers. At the same total Lya absorption 
strength some clouds have strong metal lines while others show 
no trace of any metal absorption. Further studies of larger sam- 
ples will be needed to determine the physical cause of these 
differences. 
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